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Abstract. Application of thermodynamics to driven systems is discussed. As particular examples, simple 
traffic flow models are considered. On a microscopic level, traffic flow is described by Bando's optimal 
velocity model in terms of accelerating and decelerating forces. It allows to introduce kinetic, potential, 
as well as total energy, which is the internal energy of the car system in view of thermodynamics. The 
latter is not conserved, although it has certain value in any of two possible stationary states corresponding 
either to fixed point or to limit cycle in the space of headways and velocities. On a mesoscopic level of 
description, the size n of car cluster is considered as a stochastic variable in master equation. Here n = 
corresponds to the fixed-point solution of the microscopic model, whereas the limit cycle is represented by 
coexistence of a car cluster with n > and free flow phase. The detailed balance holds in a stationary state 
just like in equilibrium liquid-gas system. It allows to deflne free energy of the car system and chemical 
potentials of the coexisting phases, as well as a relaxation to a local or global free energy minimum. In this 
sense the behaviour of traffic flow can be described by equilibrium thermodynamics. We flnd, however, that 
the chemical potential of the cluster phase of traffic flow depends on an outer parameter - the density of 
cars in the free-flow phase. It allows to distinguish between the traffic flow as a driven system and purely 
equilibrium systems. 

PACS. 05.70.-a Thermodynamics - 05.10.Gg Stochastic analysis methods - 89.40.-a Transportation 



. 1 Introduction 

, An extension of thermodynamic concepts from equilib- 
■ rium to nonequilibrium or driven systems is one of the 
' fundamental problems in physics. It refers also to so- 
called nonphysical systems, like traffic or granular flow, 
economics, biological systems, etc., where the laws of mi- 
croscopic interaction and motion differ from those known 
' in physics. Different approaches have been developed till 
I now. In the geometrical formulation of thermodynamics [Ij , 
. the latter is regarded as a theory arising in the analysis of 
' dynamics. In this concept the equilibrium thermodynam- 
ics is represented by a manifold of time-independent equi- 
'■ librium states, whereas the thermodynamics of driven sys- 
tem is represented by a manifold of slowly evolving states. 
A fc-component system undergoing chemical reaction is 
considered as an example in [1] . A more widely discussed 
approach is based on the introduction of entropy and 
usage of the entropy maximization principle in various ap- 
plications, e. g., linear dissipative driven systems [3] and 
single-lane traffic [2J • An appropriate definition of temper- 
ature is a relevant question when we speak about a non- 
physical system. In [2j the temperature T and pressure p 
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of traffic flow have been introduced via derivatives of cer- 
tain thermodynamic functions, and it has been found that 
T is negative at typical velocities. In another approach [4] 
similarities between traffic and granular flow have been 
discussed proposing two effective temperatures: one char- 
acterising fast or single-car dynamics, and another — slow 
or collective dynamics of traffic flow. 

As mentioned in [2], entropy need not occupy a posi- 
tion of primacy in a general theory beyond the classical 
equilibrium thermodynamics. We have found that in cases 
where the stationary state of a driven system has the prop- 
erty of detailed balance in a space of suitable stochastic 
variable, the thermodynamic potential can be easily intro- 
duced based on this property in a complete analogy with 
equilibrium systems. This approach can prove to be useful 
in many applications due to its relative simplicity. As an 
example we consider formation of a car cluster in one-lane 
traffic and show its analogy with the phase separation in 
supersaturated vapour-liquid system. 

The aggregation of particles out of an initially homoge- 
neous situation is well known in physics, as well as in other 
branches of natural sciences and engineering. The forma- 
tion of bound states as an aggregation process is related 
to self-organization phenomena [51[5117] . The formation of 
car clusters (jams) at overcritical densities in traffic fiow 
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is an analogous phenomenon in the sense that cars can 
be considered as interacting particles [SllQlfTO], The devel- 
opment of traffic jams in vehicular flow is an everyday 
example of the occurrence of nucleation and aggregation 
in a system of many point-like cars. For previous work fo- 
cusing on the description of jam formation as a nucleation 
process, see [TT | ll2 l [T3 j [T4]. It is related to phase separa- 
tion and metastability in low-dimensional driven systems, 
a topic which has attracted much recent interest [TSlfTBl 
[T7 l fT8 j fT9] . Metastability and hysteresis effects have been 
observed in real traffic, see, e. g., 20,21,22,23,24,25,26] 
for discussion of empirical data and the various different 
modelling approaches. 

Here we focus on the application of thermodynamics 
to such a many-particle system as traffic flow. In a first 
step we do not consider real traffic with its very com- 
plicated behaviour but limit our investigations to simple 
models of a directional one-lane vehicular flow. We hope 
this will trigger further development to describe more re- 
alistic situations of multi-lane traffic as well as of synchro- 
nized flow ^24j . We have found certain analogy with phys- 
ical systems like supersaturated vapour-liquid, although 
there are also essential differences, since the traffic flow is 
a driven system. We would like to outline some basic ideas 
and concepts developed throughout the paper. 

1. On a microscopic level traffic flow can be described 
by Bando's optimal velocity model. In this case the 
equations of motion can be written as Newton's law 
with accelerating and decelerating forces and one can 
define the potential V and the kinetic T energy of the 
car system, as well as the total energy E — T + V . 
The latter one has a thermodynamic interpretation as 
(E) = U, where U is the internal energy of the system. 

2. Traffic flow is a dissipative system of driven or active 
particles. It means that the total energy is not con- 
served, but we have an energy balance equation 



dE 
It 



<P = 



with the energy flux <P following from the equations of 
motion and consisting of dissipation (due to friction) 
and energy input (due to burning of petrol). 

3. In the long-time limit the many-car system tends to 
certain stationary state. In the microscopic descrip- 
tion it is either the fixed-point or the limit cycle in 
the phase space of velocities and headways depending 
on the overall car density and control parameters. The 
stationary state is characterised by certain internal en- 
ergy. 

4. On a mesoscopic level traffic flow can be described by 
stochastic master equation, where stochastic variable 
is the number of congested cars n, i. e., the size of 
car cluster. In this case the fixed-point solution corre- 
sponds to n = 0, and the limit cycle — to coexistence 
of a car cluster with n > and free flow phase. 

5. In the space of cluster size, the detailed balance holds 
for the stationary solution just like in equilibrium phys- 
ical systems. It allows to describe various properties of 
the stationary state by equilibrium thermodynamics. 



In particular, we calculate free energy of the system 
and chemical potentials of coexisting phases in a com- 
plete analogy with the known treatment for a super- 
saturated liquid-gas system. 
6. In distinction to equilibrium systems, the chemical po- 
tential of the cluster phase of traffic flow is not an in- 
ternal property of this phase, since it depends on an 
outer parameter - the density of cars in the free-flow 
phase. It allows to distinguish between the traffic flow 
as a driven system and purely equilibrium systems. 



2 Microscopic optimal velocity model of 
traffic flow 

Traffic flow can be viewed as a random dynamical sys- 
tem [271128] of active or intelligent particles t29|..30l[3T] . To 
describe it on a microscopic level, here we use a simple ver- 
sion of Bando's optimal velocity (OV) model for point-like 
cars moving on a one-lane road with periodic boundary 
conditions. The model is defined by the following set of 
equations 



dvj 
'dt 

dxi 

It 



1 



(Wopt(Z\Xi) - Vi 



(1) 

(2) 



where the coordinate Xi [t) as well as the velocity Vi (t) of 
each car i = 1, . . . , iV at every time moment t can be cal- 
culated out of the initial values by integrating the coupled 
equations of motion. Here 



{Axy 



(3) 



is the optimal velocity function depending on the headway 
distance Axi = Xi+i — Xi proposed in ^11. 12. 13,35, .36]. It 
includes the maximal velocity Vmax and the interaction 
distance D as parameters. Eq. ^ can be written as 



where 



m 



,(ui) + Fdec{Axi) 



(4) 



Facc{vi) = — (Vmax - Vi) > (5) 
T 
7TL 

Fdec{AXi) = — (VoptiAxi) - Vjnax) < (6) 



are the accelerating and decelerating forces, respectively. 
Similar representation has been introduced already in [9l 
ilOj . The only distinguishing feature is that in our case 
the deceleration force is specified by whereas in 
it is related to a power-like interaction potential. The 
coordinate-dependent force term is due to interaction be- 
tween cars 



FdeciAx) 



(Ax) 



T \D^ + {Axy 



- 1 
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and is always negative, starting at Fded^x = 0) = 
—Vmax'm/T, approacliing zero at infinite distances. The 
potential energy of car system can be defined as ^ = 
"Ylu^i 4'{^Xi), where 0(Z\xi) is the interaction potential of 
the i-th car with the car ahead, which is given by 



d(l){xi+i - Xi) _ d(l){Axi) 



dxi 



dAxi 



By integrating this equation we get 



(t){Ax) 



Dm 



TV 1 




arctan 




2 \ 





(8) 



(9) 



where the integration constant is chosen such that 
4>{oo) = 0. For comparison, the interaction potential of 
the form (j){Ax) oc (zix)~" has been considered in [9]. Note 
that Fdec{Axi) is not given by —dV/dxi, since the latter 
quantity includes an additional term —d4>{xi — Xi^i) / dxi. 
This term is absent in our definition of the force because 
the car behind does not influence the motion of the actual 
i-th vehicle. It reflects the fact that, unlike in physical 
systems, the third Newton's law does not hold here. 
The total time derivative of the potential energy is 



dV_ 
'dt 



N r 



E 



1=1 

N 



d(f){Axi) dxi d4>{Axi) dx 



dxi dt dxi+i 

^{Vi+i - Vi)Fdec{AXi) 



dt 



(10) 



The total time derivative of the kinetic energy 
T — is obtained by multiplying both sides 

of ^ by Vi and summing over i. It leads to the following 
energy balance equation 



(11) 



for the total energy E ^ T + V oi the car system, where 



N 

E 



[ViFacciVi) + Vi+iFdec{AXt)] 



(12) 



is the energy flux. It includes both energy dissipation due 
to friction and energy input from the engine. Eq. (|11|) 
shows that, in distinction to closed mechanical systems, 
the total energy is not conserved in traffic flow. Never- 
theless, it approaches a constant value in the long-time 
limit, where the system converges to one of two possible 
stationary states: either to the fixed point Axi = Axhom, 
Vi = Vopt (Axhom) (where Axhom = L/N is the distance 
between homogeneously distributed TV cars over the road 
of length L), or to the limit cycle in the phase space of 
headways and velocities. Both situations are illustrated in 
Fig. [TJ At small enough density of cars there is a stable 
fixed point (solid circle), which lies on the optimal velocity 
curve (dotted line). An unstable fixed point (empty cir- 
cle) exists at larger densities. In the latter case any small 
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8 
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Fig. 1. Fixed points (circles) and limit cycle (solid line) in 
the space of headways Ax and velocities v of cars. The solid 
circle represents the stable fixed point at the car density p = 
N/L = 0.0303 m~^. The empty circle is the unstable fixed 
point at a larger density p = 0.0606 m~^, where the long-time 
trajectory for any car is the limit cycle shown. The fixed points 
lie on the optimal velocity curve (dotted line) given by (|3} . The 
parameters are chosen as = 60, D = 33 m, Vmax ~ 20 m/s, 
r = 1.5 s, and m = 1000 kg. 
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Fig. 2. The total energy E of the car system, measured in units 
of mv'^ax/'^, depending on time t given in seconds. The same 
sets of parameters have been used as in Fig. [1] The upper solid 
line corresponds to a larger density p — 0.0606 m~^ where the 
limit cycle forms, whereas the lower dashed line — to a smaller 
density p = 0.0303 m~^ where the convergence to stable fixed 
point is observed. 



perturbation of the initially homogeneous fixed point sit- 
uation leads to the limit cycle (solid line) in the long-time 
limit. 

In the thermodynamic interpretation the mean energy 
(E) is the internal energy U of the system. The latter one 
thus has certain value in any one of the stationary states. 
The temporal behaviour of E for the same sets of parame- 
ters as in Fig.[T]is shown in Fig. [21 In the case of the conver- 
gence to the limit cycle (solid line) for p = 0.0606 m~^, one 
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can distinguish 6 plateau in the energy curve. The first one 
represents the short-time behaviour when starting from 
an almost homogeneous initial condition with zero veloc- 
ities, and the second plateau is the unstable fixed point 
situation. Further on, 4 car clusters have been formed in 
the actual simulation, and this temporal situation is rep- 
resented by the third relatively small plateau. The next 
three plateau with 3,2, and finally 1 car clusters reflect the 
coarse graining or Ostwald ripening process. The dashed 
line shows the convergence to the stable fixed point value 
at p = 0.0303 m-i. 

Apart from the internal energy, other thermodynamic 
functions can be introduced as well. In the following sec- 
tions we will calculate the free energy F of the traffic flow. 
By using the known relation F = U — T* S we can calcu- 
late also the entropy S of traffic flow for a properly defined 
'temperature' T*. 

Up to now we have considered purely deterministic 
equations of motion. Randomness can be included, e. g., 
by adding a multiplicative noise term to (j4|). It leads to 
stochastic differential equations 

mdv^it) = {Facc{v^) + Fdec{^X,))dt + aV^dW^{t),{13) 

dxi{t) = Vidt (14) 

considered and solved numerically in [13| . Here a is the 
noise amplitude, and dWi{t) is the increment of Wiener 
process. Similar equation with additive noise term has 
been studied in [QlllOj. An advantage of the version with 
multiplicative noise is that it guarantees the positiveness 
of velocities Vi. In the deterministic model the departure 
(leaving a cluster) times are strongly correlated in such 
a way that, in the stationary regime, one car leaves the 
cluster after each time interval of a given length ti. The 
arrival (adding to a cluster) times also are strongly corre- 
lated due to the repulsive forces. The noise makes these 
correlations weaker. It allows to apply the formalism of 
stochastic Markov processes to describe approximately the 
fluctuations of the cluster size, as discussed in the follow- 
ing section. 



3 Mesoscopic stochastic model of traffic flow 

It is easier to study the formation of a car congestion on a 
mesoscopic level, as it has been done in |ll f.l2 . 13.35) j36]. 
where we do not follow each individual car, but only look 
for the number of congested cars n, i. e., the size of car 
cluster. In this description it is also very easy to introduce 
the randomness, by considering n as a stochastic variable. 
Following [TTlfT2l[T3[l35il36] , in the simplest model only one 
cluster on a circular road is considered, and the probability 
p(n, t) that it contains n cars at time t is given by one-step 
master equation 



dpin, t) 
Jt 



— w+in — l)p(n — 1, t) 

+ w^{n + l)p(n + l,t) 

— [w+{n) + 7«_(n)] p{n,t) 



(15) 



For n — and n — N the equations look different, i. e., 
terms with p(— 1, t) and p{N + 1, t) are absent. In the sim- 
plest version of the model of point-like cars the transition 
rates are ly- (n) = 1/r and it;+(n) = Vopt{Axfree) /^a^/ree, 
where r is a reaction time constant and Axfreein) = 
L/{N — n) is the mean headway distance in the free flow 
phase. In this model no large stable cluster forms at low 
densities of cars, whereas a macroscopic fraction of them 
are condensed (jammed) into the cluster above certain 
critical density. The first situation corresponds to the fixed- 
point solution of the Bando model, whereas the second 
one — to the limit cycle. It is a remarkable fact that the 
stationary solution p"* (n) = \\mt~^oaP{n,t) obeys the de- 
tailed balance condition p^*-{n) w+ (n) = p^* (n -t- 1 ) w_ (n + 
1). It allows to describe some properties of the model 
by equilibrium thermodynamics in analogy to the liquid- 
vapour system, as discussed in the following sections, in 
spite of the fact that the traffic flow is a driven, i. e., 
nonequilibrium system. 



4 Free energy of the liquid-gas system 

The principle of detailed balance is useful to describe the 
equilibrium in a physical system. In particular, we analyse 
the condensation of supersaturated vapour to show how 
the free energy and chemical potentials can be derived 
based on this principle with the following idea to apply 
the same scheme for traffic flow. 

For simplicity, we consider a situation where only one 
cluster of molecules coexists with the vapour phase. The 
number n of molecules called monomers binded in the clus- 
ter is a stochastic variable, whereas their total number 

in a given volume V is fixed. The stochastic events 
of adding or removing one monomer are characterised by 
transition rates w+ (n) and W- (n) depending on the ac- 
tual cluster size n. Following [7[|13|. the detailed balance 
reads 



'+in-l ) 
i/7_ (n) 



= exp 



F{n) - F{n - 1) 



(16) 



where T is temperature, ks the Boltzmann constant, and 
F{n) is the free energy of state (including all possible mi- 
croscopic distributions of coordinates and momenta of free 
monomers) with cluster size n. For large enough n (|16p can 
be approximated as 



w+{n) 



exp 



W- (n) 

which leads to the equation 



dF/dn 



In 



From this we get 



W-{n) 



1 dF 



ksT dn 



F^Fn- kuT 



<n< N 



1 In 


'w+{n'y 


w_ (n') 





dn' , 



(17) 



(18) 



(19) 
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where Fq = F{n = 0) does not depend on the cluster size 
n. It is the free energy of the system without cluster, in 
this case the free energy of an ideal gas. We insert here 
the physical ansatz for the transition rates (see [13] ) 



w+{n) ^ Ag(T)(jV- n) 

W-{n) V 



■ exp 



fn-l{T)-UT) 



(20) 

where V is the fixed volume of the system, fn{T) is the 
binding energy of a cluster of size n, 
Ao(T) = /i/(27rmfcBT)i/2 is the de Broglie wave length 
of a monomer, and h is the Planck's constant. By using 
the approximation fn-i{T) — fn{T) — — 9/„(T)/9n, we 
obtain 



ksT J In 





^\UT)iN~n') 



V 



dn' + U{T) . (21) 



The integration, using ^ \nx dx = xlnx ~ x, yields 



F ^ Fo~ ksTN 
+ kBT{N~n) 



In (^Ag(T)^ 
In (xliT) 



N 
V 

N ~n 



(22) 



V 



fn{T) 



The free energy of ideal system (gas) Fq cannot be ob- 
tained from the detailed balance relation. It is given by 
Ff) = —ksThiZid, where Zid is the partition function of 
the ideal gas 



3Af 



Zid — 



-n- 

a=l 



dxa / dpa exp 



Po 



2mkBT 



1 



N\ \XoiT)J 



-oo 
3N 



(23) 



Hence, applying the Stirling formula In A^! ~ TV In — A^, 
we obtain 



Fo = keTN 

By inserting (|24p into 
sion 



ln(Ag(T)^)-l 



(24) 



F = kBT{N - n) 



In Ai^(T) 



we recover the known expres- 

N -n^ 



V 



- 1 



MT) (25) 



for the free energy of liquid~gas system under isothermal 
and isochoric conditions. The binding energy fn{T) can 
be written as 



fn{T) = ^.oa{T)n + aA{n) 



(26) 



where ^,oo{T)n represents the volume contribution, 
l^oo{T) < being the chemical potential for a flat droplet 
interface (with infinite radius r), and aA{n) is the surface 
contribution. Here cr > is the surface tension, whereas 



A{n) — Airr^ is the surface area of a droplet with ra- 
dius r. Taking into account that the number of particles 
(molecules) in the cluster is n = (ccjust 47r/3) r'^, where 
Cciust is the particle density inside the cluster, Eq. (|22p 
can be written as 



1 



n 
N 



F-Fo ff 

vk^ = n[ 

3 

+ 2 



In 1 



n 

keT N 



+ 1 



(27) 



2/3 



where p = N/V is the overall density and 



CclustkBT 



(28) 



is the diffusion length (width) of the liquid-gas interface. 

Further on we introduce dimensionless density p — 
Xq{T)p and dimensionless volume V = V/Xq{T). In this 
notation the equation ([^0)1 transforms to 



w+{n) / n 
W-[n) V A* 



PoojT) 

ksT 



(29) 



xexp 



[-m{c.u.^^iif''v-^'-'p-'''{^)"''') 



whereas (P7| becomes 



VkeT 



= P 



Poc{T) n 



ln(l--l-l 



keT N 

1/3 t7-1/3~-1 



(30) 



2/3 



These equations allow us to calculate the ratio 
w+(n)/?i'_(n), as well as the normalised (dimensionless) 
free energy difference {F — Fq)/ {VksT) depending on the 
fraction of condensed molecules n/N at a given overall 
density for fixed volume and temperature. The results of 
calculation for three different dimensionless densities p — 
5 • 10"^, 10~^, 1.2 • 10~^ at the values of dimensionless con- 
trol parameters P-oo/{kBT) — —12 and 
^(T)(cci„st47r/3)^/^F-i/3 = 0.003 are shown in Figs. [3] 
and m Note that the extrema of F — Fo in Fig. |4] cor- 
respond to the crossing points with the horizontal line 
wj^{n)/w-{n) = 1 in Fig.[3l At the smallest density (dot- 
dashed line) there are no crossing points and the free en- 
ergy is a monotonously increasing function of n/N ^ show- 
ing that the stable state of the liquid-gas system contains 
no liquid droplet. Stable droplet appears at larger densi- 
ties (dashed and solid lines) by overcoming a nucleation 
barrier (local free energy maximum in Fig. [4]) . 

The parameters we have chosen arc quite realistic, i. e., 
comparable with those of water at T = 300 K and V — 
5 ■ 10~^^ m^ with about 37 250 molecules (mass m = 2.99 • 
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Fig. 3. The ratio of transition rates ui+(n)/u)_(n) depending 
on the fraction of condensed particles n/N for three dimension- 
less densities j5 — 5 ■ 10~^ (dot-dashed line), p = 10~^ (dashed 
line), and p = 1.2 ■ 10"^ (solid line). 




0.2 0.4 0.6 jyN 



Fig. 4. Normalised free energy difference (F — Fo) / (VkBT) = 
f — fo depending on the fraction of condensed particles n/N 
for three dimensionless densities p — 5-10"'^ (dot-dashed line), 
p = 10^^ (dashed line), and p = 1.2 ■ 10~^ (solid line). 



10-23 kg) at p = lO-'^. For water at T = 300 K we have 
Ao(r) = 2.377 • 10-" m and Cdust = 3.346 • lO^^ m-^. 
Hence, the dimensionless density in the cluster pciust = 
Cciust^o — 4.491 • 10-* exceeds about 50 times the critical 
mean density p ~ pc — 9.2-10"^ at which the condensation 
(i. e., a minimum of free energy at n/N > 0) appears in our 
calculation. Assuming the above parameters of water, we 
obtain i{T) = 8.953 • IQ-^'' m for the width of the liquid- 
gas interface and surface tension a = 6.20 • 10® Nm/m^. It 
is about 3 times the characteristic intermolecular distance 
in the cluster, which is roughly c^;^{j — 3.1 • 10^^*^ m. 
The critical density pc increases with temperature and 
becomes equal to the cluster density pdust at the critical 
temperature T — T^. In our description the physically 
meaningful densities are restricted by p < pdust ■ It means 
that no condensation phase transition takes place for these 
physical densities at T > Tc. Assuming that /ioo and a do 
not change with temperature and that the above given 
values of dimensionless control parameters correspond to 
T = 300 K, we find ~ 430 K in our example. 



suming the model of point-like cars we have Ax free (n) = 
L/{N — n), where L is the length of the road and N is the 
total number of cars. It yields 



(n) 



1 - n/N 



1 + p^D^{l-n/Ny 



(32) 



where p =^ N/L is the car density. Introducing the di- 
mensionless density p = pD and a dimensionless control 
parameter b = D / (vmaxT) , it becomes 



In 



w+{n) 
W- (n) 

w+{n) 



1 p{l - n/N) 
b l + p'^{\-n/Nf 



(33) 



w_(n) 



In 



In 



In 



n 
N 



-2 f. n 
P V-N 



(34) 



5 Free energy of traffic flow 

Now wc make similar calculation of free energy for the 
traffic flow model introduced in Sec. [31 Similar general re- 
lations (fT6)) to p9)) are valid to describe the stationary 
(quasi-equilibrium) properties in the space of car cluster 
size n, since the detailed balance is the property of the 
stationary solution of the one-dimensional one-step mas- 
ter equation for the probability distribution over n. Here 
we only replace /cgT with T* which is the 'temperature' 
of traffic flow having energy dimension. 

The ratio of transition rates in this case reads 



w+ [n) 
W-{n) 



fopt 



{Axfr 



Aa 



free 



(31) 



where Vopt (Ax) is the optimal velocity function given by ^ 
and Ax free is the headway distance in free flow phase. As- 



By inserting the latter relation into (fTO)) (where fc^T 
T*), the integration using / In (l + a;^) = 2 arctanx -I- 



a; In (1 + x^) — 2x yields 

F-Fo 



LT* 



P 



1 



n \ , f n 



n 
N 



n 



~2 
P 



n 



1 

N 



In (1 



^ 1 



P'] 



-2 arctan p — 2 arctan ( p 



71 ■ 
1 

N. 



(35) 



where L = L/D is the dimensionless length of the road. 

The results for w+{n)/w-{n) and {F-Fq)/(LT*) de- 
pending on the fraction of congested cars n/N at four 
different densities are shown in Figs. [5] and [SI The value 
of the dimensionless control parameter has been chosen 
& = 2/7 « 0.2857. It corresponds, e. g., to £> = 24 m, 
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n/N 



Fig. 5. The ratio of transition rates w+{n)/w-{n) depending 
on the fraction of congested cars n/N for four dimensionless 
densities p = 0.1 (dotted hne), p = 1 (dot-dashed line), p — 
3.186 (dashed line), and p = 5 (solid line). 




n/N 



Fig. 6. Normalised free energy difference {F — Fo) / (LT*) — 
f ~ fo depending on the fraction of congested cars n/N for 
four dimensionless densities p = 0.1 (dotted line), p = 1 (dot- 
dashed line), p = 3.186 (dashed line), and p = 5 (solid line). 



only shifted mean value (p) = m{v) = mvopt {'^Xhom) in 
accordance with the optimal velocity Vopt i^Xhom) in the 
homogeneous flow of cars with the mean headway distance 
^Xhom = L/N = 1/ p. The Gaussian form of the distri- 
bution is well consistent with the simulation results for 
the stochastic car-following models [9l fT0[|T3] . We should 
take into account also that cars are moving always in one 
direction, i. e., momentum p > Q always holds. Finally, 
in distinction to the ideal gas of molecules, the coordi- 
nates and momenta of cars are one-dimensional. Hence, 
by analogy to we can write 



N 



1 / L 



2mT* 





N 



m \XoiT*) 



(36) 



where Ao(T*) = h/{2'KmT*y/^. The latter approximate 
equality in holds when (p)'^ / {2mT*) ^ 1 or, in other 
words, when the width of the velocity distribution is nar- 
row as compared to the mean velocity. The latter con- 
dition is well satisfied for the model (fT51) — ([T^ with 
certain set of control parameters used in the simulations 
of [13] (see Fig. 40 there). The distribution width, how- 
ever, increases with the noise amplitude. In fact, the ap- 
proximation (j36[) is good enough when the distribution 
function has small value at zero momentum p = 0, as in 
the simulation results of [51llO|. According to the above 



consideration, temperature in traffic flow is a parameter 
which controls this distribution width or the amplitude 
of velocity and momentum fluctuations, like in the ideal 
gas of molecules. According to (|36t . the ideal part of free 
energy reads 



Fo = -T* [N\n{L/Xo{T*)) - \nNl] 
~r*iV[ln(pAo(T*))-l] . 

This expression is analogous to (|24p . 



(37) 



Vmax ~ 42 m/s, and r = 2 s. Like in the case of the liquid- 
supersaturated vapour system, w^{n)/w-{n) is never 1 
and no stable car cluster forms at small densities (dotted 
line). In distinction to the liquid-vapour case, the cluster 
appears without a nucleation barrier in the actual traf- 
fic flow model at somewhat larger densities p (dot-dashed 
line), whereas the nucleation barrier (free energy maxi- 
mum) shows up only at even larger p values (solid line). 

In the above calculation we have determined only the 
difference F — Fq, but not the free energy Fq of the ideal 
system without car cluster. Like in the case of supersat- 
urated vapour, the latter cannot be derived from the de- 
tailed balance. It should be calculated from a microscopic 
model. Now, however, we should take into account that 
the distribution over momenta for cars is not the same 
as that for molecules in an ideal gas. As a first approxi- 
mation we may assume similar Gaussian distribution with 



6 Relaxation to a free energy minimum 

Now we consider the general behaviour of a system (equally 
valid for liquid-gas system with T* = kgT and traffic 
flow) in vicinity of a local or global minimum of F{n). In 
this case the argument of exponent in (|17p is small and 
we can make a Taylor expansion 

^.expf-^Ul-^^. (38) 
w-{n) ^\T*J T* dn ^ ' 



It can be rewritten as 



w^{n) dF 



(39) 



On the other hand, we can write in a deterministic ap- 
proximation 

— = w+{n) — w-{n) . (40) 
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Comparing ([39| and (|40|) . we obtain 



dn 



w-{n) dF 
T* 9^ 



(41) 



Like in the Landau theory of phase transitions, we can 
expand the free energy around the minimum point n = no 
defined by 

OF 

TT = ■ (42) 

In the first approximation, where we retain only the lead- 
ing term, we have also w+(ri) = w-{n) = w± (no). It leads 
to the kinetic equation 



dn 
It 



where 



To 



~ -Fq (n - no) , 
w±{no) d^F 



d'n? 



(43) 



(44) 



is the relaxation rate. For F^ > 0, what corresponds to 
minimum of F, the solution is the exponential relaxation 
to 71 = no, i. e.. 



n(t) = no + (n(0) — no) e 



-To t 



(45) 



This solution is valid also for Fq < 0, in which case no 
corresponds to a free energy maximum. In this case it 
describes the deviation from this maximum point. 



7 Chemical potentials 

Our system can be considered as consisting of two phases: 
the cluster phase with n particles and free energy Fci{n), 
and the ideal gas phase with Ni^ — N — n particles and 
free energy Fid [Nid)- The total free energy then is F = 
Fci + Fid- While the total number of particles TV is fixed, 
the number of particles in any of the phases fluctuates. Ac- 
cording to the definition, we can write fXci = dFd /dn and 
Hid = dFid/dNid = —dFid/dn for the chemical potentials 
of these phases. Hence 



dF 
dn 



OF,, dF, 



id 



dn 



dn 



and the kinetic equation (|4ip can be written as 



dn w± (no) 
It ~ T* 



(46) 



(47) 



The latter equation has certain physical interpretation: 
the driving force pushing the system to the phase equi- 
librium is the difference of chemical potentials in both 
phases. The equilibrium is reached when the chemical po- 
tentials of the coexisting phases are equal, i. e., fid = tJ-id- 
For the hquid-gas system Fid{T,V,N,n) is given by 
where N is replaced with Nid = N — n, i. e., 



F,d{T,V,N,n)^kBT{N -n) 



In Ai*(r) 



N- 



V 



1 



(48) 



Hence the total free energy (l^5|) can be written as 

F = F,d{T,V,N,n) + fniT). (49) 

The chemical potential of the liquid phase is thus given by 
the derivative of the binding energy /n(T) = Fci{T, V, N, n), 
i. e., 

^lcl ^ t^ooiT) + a^^ = fi^{T) + kBTe{T)k{n) , (50) 



dn 



where 



k{n) = - ^ {cciust 47r/3)'/'n-i/3 



(51) 



is the curvature of the liquid surface for a droplet with 
radius r and surface area A(n) = Airr"^ . The chemical po- 
tential of the gaseous phase calculated from is 



Mid = 



^^"=fcBTlnfA3(T)^-" 



dn 

kBT\n{pgas) , 



V 



(52) 



where pgas — \^){T){N ~n) /V is the dimensionless density 
of molecules in the gaseous phase. According to these ex- 
pressions for the chemical potentials, the ansatz ((20|) can 
be written as 



w+ (n) 
w_(n) 



( ^J^^d\ f fn-l{T) - fn{T) 



ksT 



exp 



keT 



(53) 



The latter relation is consistent with (fT7|) and (|46|) . 

By analogy, the free energy of the free flow phase in 
traffic is 

" / N - n\ 

Fu{T\L,N,n)^T*{N -n) hi f Ao(r*)— -1 

(54) 

as consistent with ([37|) where we put N Nid — N — n 
and p — > Nid/L. From ([54]) we get 



(55) 



where pfree = D{N — n)/L is the dimensionless density 
of cars in the free flow phase. The chemical potential of 
the cluster phase can be easily calculated from Eqs. ([18]), 
([341) . and It yields 



Mc; = ~T* { In 
= -T* <! In 



D 

D 
A 



-In 



1 + 1 _ _ 



-ln[l + pj_]| 



(56) 



It is remarkable that, in distinction to the liquid-gas sys- 
tem, the chemical potential of the cluster phase is not 
an internal property of this phase, since it depends on the 
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outer parameter - the density of the surrounding free-flow 
phase p free ■ The physical interpretation of this fact is that 
the traffic flow is a driven system, which approaches a sta- 
tionary rather than equilibrium state in the usual sense. 
However, as we have shown here, various properties of this 
stationary state can be described by equilibrium thermo- 
dynamics. 

Free energy of the cluster phase can be calculated con- 
sistently from , (|57)l , and ([51]) according to F = Fei -\- 
Fid- The result is 



Fei{T*,L, N, n) = T*N ^ (2 + In 



arctan p — arctan ( p 



1 - 



n 

1 

N. 

n 



(57) 



N 



In (1 + p2 



It is consistent with /ic; ~ dFei/dn. 



8 Conclusions 

1. In the current paper we have shown how thermody- 
namics can be applied to such a many-particle system 
as traffic flow, based on a microscopic (car-following) 
as well as a mesoscopic (stochastic cluster) descrip- 
tion, in analogy to equilibrium physical systems like 
supersaturated vapour forming liquid droplets. The 
basic idea here is to derive the free energy function 
and chemical potentials by using the detailed balance, 
which holds in the stationary state of traffic flow in the 
space of car cluster sizes. 

2. Distinguishing features between the traffic flow and 
equilibrium physical systems have been outlined. In 
particular, we have found that the third Newton's law 
does not hold on the level of " microscopic" equations 
of motion for individual cars. Besides, the traffic flow 
is a dissipative system with inflow and outflow of to- 
tal energy. Unlike in equilibrium systems, the chemical 
potential of the phase of congested cars is not an in- 
ternal property of this phase, since the traffic flow is a 
driven system. 
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